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Abstract 

This paper presents a new method for the regularization 
of the inverse problem of electrocardiography. Our goal is 
to reconstruct the epicardial potentials (EPs) from mea- 
sured body snrfaoe potentials. Tliis non-invasive tedi- 
nique is usefiil for the diagnostic of the cardio-vascular 
diseases. In our approach, the time-correlation between 
EPs is \ised as the regularizing a jmori information. This 
information is introduced via the state-space representa- 
tion proposed in [I]. Our contribution is the derivation of 
maximum likelihood estimators for identification of the 
parameters of the state-space model, and in a second 
stage, fior determination of the EPs. Therefore, all un- 
known quantities are determined &om the only available 
data: the body aur&ce potentials. 

1 Problem statement 

The goal of the inverse problem of electrocardiography 
ia to quantitatively evaluate the cardiac electrical activity 
from measured torso potentials. The relationship between 
body surface potentials and EPs that can be considered 
linear with a good approximation. Therefore, at time- 
sample n, one can write 

iv{n)=H^B{n) + u{p) n^l...N, (1) 

where and respectively denote the vectors of body 
surface potentials and EPa, and where is the trans- 
fer matrix between the heart and the torso; w is a sero- 
mean Gaussian white noise vector with covariance matrix 
rl{l is the identity matrix) which represents modeling 
and measurement errors. In this study, H was computed 
using a 3-D finite element model of the torso according to 
the technique described in [2]. 

It is well-known that this problem is ill-posed, and that 
it must be regularized through incorporation of prior in- 
formation on the solution. In the clasdcai regulariza^ 
tion methods (Tikhonov regularization or regularization- 
truncation) [3], the time-correlation of the cardiac activa- 
tion process is not accounted for as the data are indepen- 
dently processed time-frame by time-frame. 


Oster & Rudy [4] pointed out the importance of tem- 
poral information for regularization, but they did not ex- 
plain how to account for this information in a non-invasive 
manner. More recently, Brooks & al. [5] proposed a par- 
tial answer, but their fbrmnlation of the a priori informa- 
tion lacks flexibility and the estimation of the regulariza- 
tion parameters remains a problem. 

In order to account for the time-dependence across 
tfane-frames, Jdy te al. [1] proposed the foUowing time- 
invariant linear prediction equation: 

*B(n) = r#E(n-l)-l-v(n) n = l...N, (2) 

where F denotes a time-invariant transition matrix and 
V a zero-mean Gaussian white vector with covariance 
matrix qf. Experiments performeid for several patients 
showed the validity of the above model, but also indicated 
that the modd parameters vary greatly from patient to 
patient. Therefore, model parameters {F, q, r} must be 
spedfically estimated for each patient. Here, we propose 
a method for solving this problem and reconstructing the 
EPs using torso potentials only. 

2 Method 

We start with estimation of the EPs when modd (2) is 
known. With our maximum likelihood approach, the 
estimate ^s(n) maximizes Gaussian density p{^b \ 
*riJ^i9,r), or, equivalently, minimi2se8 .the following 
least-squares ctiteruin: 

N 

J = X)'|l*r(n)-Jf*,(n)!p 

-l-A||*E(n)-F*E(n-l)||», (3) 

where X = r/q. Easy computation of 4i;(n) can be car- 
ried out with a Kalman smoother. 

We now turn to maximum likelihood estimation of {F, 
g, r} from torso potentials only. The Lkelihood fiinction 
is defined by p(*r(l), ■ • ■ , *t(JV) | F, q, r). Under our 
Gaussian assumption, the likelihood can be computed re- 
cursively for any ^-alue of {F, q, r}, by usmg Kalman 
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filtering equations. However, maadmization of the lilce- 
libood presents difficulties, as no dosed-fbrm solution is 
available. Here, we propose to maxbnize the Ukdihood 
with an expectation-maximization (EM) algorithm. The 
EM algorithm is an iterative procedure that guarantees an 
increase of the likelihood at each iteration and that con- 
verges to a stationary point of the likelihood. The EM 
algorithm is well suited to our problem becatise it can 
be implemented using by-products of the same Kalmaa 
smoother as the one used for computation of #s(n) [6], 
and because it allows for various parametric forma ci ma- 
trix F. This last point a impwtant in practice because 
the number of observed data is rather limited compared 
to the number of elements in F. Therefore, parameteri- 
ution of .F is necessary so as to improve the conditioning 
of the identification problem. 

3 Results and conclusion 

EPa were measured with a 63-lead mapping system using 
a sock electrode array in 8 patients with Wolff-Parkinson- 
White s3mdrome who had undergone arrhythmia surgery. 
Each recording was taken during normal sinus rhythm and 
lasted 1.024 s. Body surface potentials were simulated 
by multiplying EPs with transfer matrix H and adding 
Gaussian noise to the results with a signal-to-noise ratio 
of20dB. 

In a first stage, several parametric forms of F were es- 
timated firom the EPs of each patient. Then, the BPs 
were reconstructed using the estimated modd and the 
simulated torso potentials using Kalman smoothing (su- 
pervised method}. Not surprisingly, we observed that EP 
reconstructions are better a^ the number parameters of F 
grows up. 

In second stage, {F, q, r} and the EPs were estimated 
for each patient firom the simulated torso potoitiab using 
the method described above. Matrix F was parameter- 
ized as: 

F = ftI + ^5. (4) 

Q and are two unknown scalars parameters; matrix 
5, known in advance, represents the time-space conela^ 
tions between EPs at any ^ven lead and their four closest 
neighb<»s at the previous time-sample. This structure re- 
lies on the assumption that propagation of the EPs is a 
local phenomeocm. 

Our results indicate that the EM algorithm yields satis- 
factory estimates of {F, r} when constraint (4) is used. 
Corresponding EP otimates e^diibit a lower reconstruc- 
tion error and a better localization of the extrema than 
TikhonoT-regularizcd solutions (see figure 1). Therefore, 
the proposed method is able to account for time-space 
correlations in an effective manner, and to improve the 
solution of the inverse problem of electrocardiography. 
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Figure 1: Reconstructed EP cards, (a) Actual EPs. 

(b) Reconstructed EPs using the supervised method. 

(c) Reconstructed EPs using Tikhonov r^ularizatiou. 

(d) Reconstructed EPs uang the proposed method. Our 
method yields better localization of the extrema and a 
lower reccmstruction error than Tildionov regularization. 
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